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Multigap superconducting state in molecular metallic hydrogen 
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The thermodynamic parameters of the superconducting state, that gets induced in the metallic 
molecular hydrogen under the influence of the pressure at 414 GPa (Tc — 84 K), have been de- 
termined. The calculations have been conducted in the framework of the three-band Eliashberg 
equations. The order parameters (A a ) and the wave function renormalization factors (Z a ) have 
been derived; the symbol a denotes the band index: a £ {a,b,c\. It has been stated that the 
dimensionless ratios 2A a (0) /fcsTc are equal to: 5.55, 3.96 and 3.53, respectively. Next, the to- 
tal normalized function of the density of states, the free energy, the thermodynamic critical field 
and the specific heat have been determined. The obtained results differ significantly from the re- 
sults achieved in the framework of the one-band Eliashberg model for the comparable value of the 
pressure (p = 428 GPa). In the last step, the values of the wave function renormalization factors 
have been estimated for Tc- It allowed us to calculate the maximum value of the electron effec- 
tive mass in a given band [([wig] ). The following results have been achieved: [m*] a = 2.99 [m e ] , 
[m* t ,] b = 2.10 [m e ] b and [m*] c = 1.94 [m e ] c , where [m e ] a stands for the electron band mass. 
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Hydrogen is an element in which the superconducting 
state with a very high value of the critical temperature 
(Tc) should be induced at the high pressure m. This 
claim is justified on the basis of the BCS expression Q, 



3] 



k B T c = 1.13wDexp - 



A 



(1) 



where fc^ is the Boltzmann constant, ujo is the Dcbye fre- 
quency, the symbol A denotes the electron-phonon cou- 
pling constant, and //* represents the Coulomb pseudopo- 
tcntial. 

In the case of hydrogen, the Debye frequency reaches 
high values due to a very small mass of the atomic nuclei 
forming the crystal lattice (single protons); the electron- 
phonon coupling is strong due to the lack of strong in- 
ternal electron shells. On the other hand, the depairing 
electron correlations modeled by the Coulomb pseudopo- 
tential are not high enough to significantly lower the crit- 
ical temperature 0], @, [I], 0- 

The studies on the thermodynamic properties of the 
superconducting state in hydrogen have lasted for many 
years. In particular, the ab initio calculations suggest 
that hydrogen becomes a subject of the metallization for 
the pressure at about 400 GPa In the range of the 



pressure from 400 GPa to 500 GPa the superconducting 
state, characterized by a high value of the critical tem- 
perature (242 K for p = 450 GPa) 0, should be induced 
in the molecular phase of the metallic hydrogen. Ad- 
ditionally, it should be noted that the superconducting 
state in the molecular phase of the metallic hydrogen can 
be strongly anisotropic (p = 414 GPa) In this case, 
an accurate description of the thermodynamic properties 
is possible only in the framework of either anisotropic 
model or multiband model [13], [ll|, fl2T |. 

Above the pressure of 500 GPa, the metallic phase 
of the molecular hydrogen becomes dissociated Q, (HI, 
fl4| . The results of the calculations performed for the 
metallic atomic hydrogen just above the dissociation 
pressure (p = 539 GPa) suggest a high critical tempera- 
ture value (Tq = 357 K); other thermodynamic parame- 
ters differ si gnifi cantly from the values predicted by the 
BCS theory [15|. 

For extremely high pressures (above 1 TPa), the super- 
conducting state properties have been analyzed in the 
works H, [13], Gil- It has been stated that for p = 2 
TPa, the critical temperature may reach a record level 
of an order 600 — 700 K. It is worth to underline that 
for p = 2 TPa, the ratio of the energy gap to the critical 
temperature (Ra) changes from 6.17 to 6.63 according to 
the assumed value of the the Coulomb pseudopotential. 
Similarly high values of the parameter i?A can be found 
only in the family of the high temperature superconduc- 
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tors 

In the presented paper, we have studied the thermo- 
dynamic properties of the anisotropic superconducting 
state, which can be induced in the metallic molecular hy- 
drogen under the pressure at 414 GPa (structure Cmca). 
In the present case, three energy gaps can be distin- 
guished, which arc connected with the different parts of 
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the Fermi surface. The first and the largest one (symbol 
a) appears on the strongly coupled disk around the T- 
point; the second gap (b) is connected with the prism-like 
sheets; the energy gap of the smallest value (c) exists on 



the remaining regions of the Fermi surface H . According 
to the above, the values of the thermodynamic parame- 
ters have been calculated with the use of the three-band 
Eliashberg equations (20j : 



/3e{a,6,c} m=-Af 



(2) 



and 



M 



Z« = l + nk B T J2 E — 



(3) 



The functions A" ee A q (iw n ) and Z° = Z Q rep- 
resent respectively: the order parameters and the wave 
function renormalization factors; u) n is the Matsubara 
frequency: u>„ = irksT (2n — 1), where fee is the Boltz- 
mann constant. The symbol a represents the electronic 
band index. In the considered case a £ {a, b, c}. 

The functions K a @ {w n — ui m ) define the pairing ker- 
nels for the electron-phonon interaction: 



(ui n - uj m )' 2 + O 2 



(4) 



In the presented work the pairing kernels have been ap- 
proximated in the agreement with the formula below: 



a/3 



(uj n —W m ) + 0% 



(5) 



where A Q/3 denote the electron-phonon coupling constants 
X afi = 2 j+oo a^(n)F ag (n) _ por the moleculax hydrogen 

under the influence of the pressure at 414 GPa, the cou- 
pling matrix takes the form [9|: 



[A' 



0.14 0.37 1.45 
0.20 0.06 0.82 
0.11 0.11 0.71 



(0) 



The symbol flc in Eq. ([5]) represents the characteristic 
phonon frequency, which value has been chosen so that 
the critical temperature estimated on the basis of the 
Eliashberg equations would be in a fair agreement with 
the value of Tc determined with the help of the ab initio 
calculations: Tc = 84 K Q. In particular, in the Eliash- 
berg equations we have assumed that T = Tc and then 
we have decreased the value of the parameter £lc until we 
have reached the equality: A° i=1 = Af„ =1 = A^ n=1 = 0. 
The resulting courses have been presented in Figure 1. 
We have found that flc equals 107.73 meV. 
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FIG. 1: The dependence of the maximum values of the order 
parameters on Qc- 



In the Eliashberg formalism, the Coulomb repulsion 
between electrons is modeled with the help of the func- 

are 



tions: fj^ (w m ) = (w c - K„|), where fj^ 

the values of the Coulomb pseudopotentials. 
responding matrix takes the form [9J: 



[ftp] 



0.046 0.010 0.100 
0.006 0.046 0.180 
0.007 0.025 0.182 



The cor- 



(7) 



The quantity u>c determines the cut-off frequency: 
ujc = Silmax, where f2 max denotes the value of the maxi- 
mum phonon frequency (O max = 496 meV). 

The Eliashberg equations have been solved for 701 
Matsubara frequency (M = 350). We have used the nu- 



merical methods presented in the works [2l| , [22( . In the 
considered case, the functions A" and Z% are stable for 
T > T = 18.9 K. 
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FIG. 2: (A)-(C) The values of the order parameters and (D)-(F) the wave function renormalization factors on the imaginary 



In Figure 2, we have presented the form of the order 
parameters and the wave function renormalization fac- 
tors on the imaginary axis for the selected values of the 
temperature. It can be easily noticed that the consid- 
ered functions have the characteristic Lorentzian shape; 
wherein the order parameters saturate much faster than 
the wave function renormalization factors. 

The temperature dependence of the functions AJ^ and 
Z^ can be investigated in the easiest way after plotting 
the courses of A« =1 (T) and Z« =1 (T) (Figure 3 (A) 
and (B)). Let us note that the curves plotted in Fig. 
3 (A) can be parametrized with the help of the function 

A" =1 (T) = A« =1 (T ) Jl^Jj^f, where the param- 
eter A™ =1 (T ) takes the values: 19.73 meV, 14.15 meV 
and 12.64 meV. The exponent ry is equal to 3.2. From 
the physical point of view, the functions 2A^ =1 (T) with 
a good approximation reproduce the temperature depen- 
dence of the energy gaps on the relevant parts of the 
Fermi surface. 

Based on the results presented in Figure 3 (B), it 
has been stated that the maximum values of the wave 
function renormalization factors very little depend on 
the temperature. Let us notice that the functions 
%m=i (T) reproduce the values of the electron effective 
mass ([ra*] n ) in a given band with a good approxima- 
tion. In particular: [m*] a ~ Z^ =1 (T) [m e ] a , where the 
symbol [j7i e ] a stands for the electronic band mass. Thus, 
the obtained results show that the effective mass of the 
electrons very little depends on the temperature in the 
whole range of the existence of the superconducting state, 



and it reaches relatively high values. 

For T = Tc, the numerical results obtained for the 
wave function renormalization factors can be compared 
with the rigorous analytical results. In particular, the 
formula below is true: 



(T C ) = 1 



Pe{a,b, 



X 



a/3 



(8) 



The numerical results have the following forms: 
Z a m=x (T c ) = 2.96, Z b m=l (T C ) = 2.08 and Z c m=1 (T c ) = 
1.93. Using the expressions © and (8), it can be proved 
that the analytical approach gives the same results, which 
demonstrates a high accuracy of the numerical analysis. 

The exact physical value of the order parameters can 
be estimated based on the knowledge of the order pa- 
rameter functions on the real axis (A Q (w)). Considered 
functions can be calculated on the basis of the formula 



A a (u) 



Pi + Pi^ 



(9) 



where p" and q" are the number coefficients; whereas 
the exponent r is equal to 7. Then, the equation below 
should be solved: 



A (T) = Re [A (w = A (T) , T)} 



(10) 



From the physical point of view, the most interesting 
result has been being obtained for the lowest considered 
temperature. Figure 4 presents the form of the order pa- 
rameters for T — Tq in the range of the frequency from 
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functions for the superconducting and normal state in the 
range of a given band. In addition, the following weights 
have been estimated: \ a = 1/3. The normalized band 
densities of states have been calculated using the formula: 



Re 



\u-tT\ 



(u> - iT) 2 - (A a (u)) 2 



(12) 



where the pair breaking parameter T equals 0.15 meV. 

In Figure 5, the form of the total normalized den- 
sity of states for the selected values of the temperature 
has been plotted. Based on the presented data, it can 
be noticed that the characteristic peaks of the function 
Ns (w) /Nn (w) appear in the points u = ±A a . 



FIG. 3: (A) Fhe maximum values of the order parameters and 
(B) the maximum values of the wave function renormalization 
factors in the range of the temperature from To to Tc- 



to 70 meV. On that basis, the values of the dimension- 
less ratios have been determined: i?A° = 2A Q (0) /ksTc, 
where A a (0) ~ A Q (To). The following results have been 
obtained: i?A» = 5.55, i?A b = 3.96 and i?A c = 3.53. Re- 
ferring to the BCS theory 0,0, which predicts that 
i?A = 3.53, one can see a particularly strong deviation 
from the BCS result for the band a. 

30 



> 

E 

v> 
t 

CD 

-I—" 

CD 

E 

Oi 

L_ 

03 
Q_ 
u_ 

CD 



20 



10 




Re[A a (co)] 

Re[A b (co)] 
Re[A c (co)] 



10 20 



30 40 

co meV 



50 60 70 



FIG. 4: The form of the functions Re [A a (w)] on the real 
axis. The functions Im [A a (cj)] in the considered area of the 
frequencies assume the value equal to zero. 

In the next step, using the results achieved for the 
order parameter on the real axis, the total normalized 
density of states has been determined: 



N N [u) 



E 



x£{a,b,c} NK ' 

where the symbols Ng (w) and Nft (w) denote the density 
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FIG. 5: The total normalized density of states in the range 
of the temperatures from To to Tc- 

After the analytical continuation of the functions Z" 
on the real axis, the exact values of the electron ef- 
fective mass for a given band can be determined. Up 
to this point, the following formula should be used: 
[m*] = Re [Z a (0)] [wie]™- On the basis of the conducted 
calculations, it has been found that the values of the ef- 
fective mass are very weakly dependent on the temper- 
ature and reach a maximum at a critical temperature. 



In particular: [m*] = 2.99 [m R |„ 



elb 



2.10 \m e ) b and 



[m*] c = 1.94 [m e ] c . The full forms of the functions Z a (u>) 
for T = Tc have been presented in Figure 6. 

The thermodynamic critical held and the specific heat 
should be calculated on the basis of the free energy dif- 
ference between the superconducting and normal state: 
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M 



AF = -2nk B T 



t£{a,6,< 



Pa (o) [y^ + tA^) 2 - M]\z^ S) z%w 



(13) 



where the upper indexes 5 and N denote the supercon- 
ducting and normal state, respectively. The vector of the 
density of states at the Fermi level takes the following 
form [911 : 
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FIG. 6: The form of Re [Z a (w)] on the real axis. In the 
considered area of the frequencies, the functions Im [Z a (w)] 
assume the value equal to zero. 

The thermodynamic critical field should be estimated 
on the basis of the expression below: 



H c = V-SttAF. 



(15) 



The specific heat difference between the superconduct- 
ing and normal state can be calculated using the formula: 



AC 



-T- 



d 2 AF 



dT 2 



(16) 



In Figure 7, we have provided the plot of the depen- 
dence of the free energy difference on the temperature. 
Additionally, the insets (A) and (B) present the form of 
the thermodynamic critical field and the specific heat dif- 
ference in the range of the temperatures from To to Tc- 
Comparing the achieved results to the results obtained 
in the framework of the one-band model for the compa- 
rable value of the pressure (p = 428 GPa), it has been 
found that the discussed functions generally reach much 
larger values. In particular for T = Tq, the following 
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FIG. 7: The free energy difference between the superconduct- 
ing and normal state as a function of the temperature. The 
inset (A) - the thermodynamic critical field as a function of 
the temperature. The inset (B) - the specific heat difference 
between the superconducting and normal state as a function 
of the temperature. 



has been achieved [A-F] p=414GPa 



6.34 



c. 



p=414GPa 



Cj p =428GPa 



/[AF] 

p=428GPa 

= 1.72; whereas 
the results are: 
4.34. From the 



and 

for the critical temperature 
[AC]p =414GPa / [AC] p=428GPa = 
physical point of view, the above results are linked with 
the relatively high value of the electron density of states 
in the bands b and c. 

Summing up, the basic thermodynamic parameters of 
the superconducting state in the metallic molecular hy- 
drogen have been provided in the presented work. The 
pressure value equal to 414 GPa has been considered. 

It has been stated that the dimcnsionless ratios 
of the order parameter to the critical temperature 
(2 A Q (0) /k B T c ) are equal to: 5.55, 3.96 and 3.53, re- 
spectively It can be easily noticed that the first value 
significantly exceeds the value predicted by the classical 
BCS model. Let us draw our attention the fact that the 
courses of the order parameters on the real axis allowed 
us to determine the total normalized function of the den- 
sity of states for the selected values of the temperature. 

Additionally, the following quantities have been cal- 
culated: the free energy difference, the thermodynamic 
critical field and the specific heat difference. On the ba- 
sis of the achieved results, it has been concluded that the 
values of the considered functions are highly underesti- 
mated in the framework of the one-band model. 
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In the last step, the wave function renormalization fac- 
tors have been determined. It ensured the calculation 
of the maximum values of the electron effective mass in 
a given band ([wi*] ). The following results have been 

] b = 2.10 [m e ] b and 



obtained: [m*] a = 2.99 [m e 
K] e = 1.94[me] c . 



port. 
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